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Synopsis 


The magnetic systems in the presence of a random magnetic field are good repre- 
sentations of many impure systems. The applicability of the model goes well be- 
yond ferromagnetic systems. The model well describes any solid state system with 
quenched impurities or lattice defects which locally break the symmetry between the 
ordered states. This is a general rule and applies to many systems in magnetism, 
ferroelectricity and suiface physics etc. [2]. 

The random magnetic fields ha\e a marked effect on the second order phase tran- 
sitions. It is now well established that the lower ciitical dimension of the random 
field Ising model is two and that of .systems with continuous symmetry is four. It is 
also clear that as far as critical behaviour is concerned, the random field fluctuations 
completely dominate the thermal fluctuations. This has led to the general consensus 
that the critical behaviour is controlled by the zeio temperature fixed point. The 
purpose of this thesis is to show the rclcvence and utility of zero temperature fixed 
point for both statics and dynamics of the random field systems. 

The thesis is divided into six chapters. The first chapter is an introduction where 
we review the existing literature. In the second chapter we introduce the phenomenol- 
ogy [16] of the zero tempeiature fixed point. It is shown that there are three indepen- 
dent exponents in contrast to two for the conventional fi.xed point. The additional 
exponent is related to the lenormnlisation of the random field or equivalently the 
renormalisation of the teniperciture. 



In the third chapter the replicat('<l approach is combined with decimation tech- 
nique to set up renormali.sation group Hows for the random field Ising model (RFIM), 
which are Migdal-KadanofF like for dimension d > 2. Considering the RFIM in 
d = 2 e, we find the critical exponents to first order in e. Our results are in 
agreement with d d — 2 + i) rule , according to which the critical behaviour of a 
d-dimensional RFIM is identical to the behaviour of pure system in d — 2 -f ?? dimen- 
sion. In d = 2, the correlation length ^ ~ e where J and Hr are the spin coupling 
and random field strength respccti\'cly, and a is an arbitrary constant. 

The system with continuous symmetry are considered in chapter four. The stan- 
dard decimation is combined with with Migdal- Kadanoff bond moving appro.ximation 
to set renormalisation flows. Considering the system in d = 4 -f- e, the critical expo- 
nents are evaluated to first order in c. The exact inequality fj < 2t] due to Schwartz 
and SofFer [15] is satisfied as equality. This would also imply that the d-dependent 
hyperscaling law du = 2 — o is modified to (d — 2 -f Tj)u = 2 — 0:. However, there is 
no exact dimensional reduction. In d = 4 the correlation length ~ e 

In chapter five we study the critical dynamics of the RFI.M. There have been 
two approaches to critical dynamics namely (i) conventional dynamics where the 
relaxation time r ~ and (ii) the activated dynamics where the relaxation time 
T ~ Recent experiments over a long range of frequencies have revealed that 
that the critical dynamics is of the activated type but the measured value of the 
dynamical exponent 0 is much smaller than the theoretical expectations. We show 
that the careful application of dynamical scaling and the zero temperature fixed point 
leads to an ever present •correction to scaling’ which causes the relevent exponent to 
differ so strongly from the theoretical expectations. 

The last chapter is devoted to one dimensional models in a random field. In section 
A, we consider the Ginzburg-Landau model in a random magnetic field. Using the 
replicated hamiltouian it is possible to reduce the problem of calculating the free 
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energy to that of calculating the the ground state energy of a quantum mechanical 
problem. The fiee energy i.s lowered and specific heat is raised beyond the respective 
zero field value. In section B, we consider the RFIM in one dimension. Considering the 
n-times replicated hamiltonian we concentrate on estimating the the largest eigenvalue 
of the corresponding 2" x 2" transfer inatri.x by using the variational principle. From 
the estimated largest eigenvalue which permits analytical continuation, we find the 
free energy for all temperatures. In section C, we take the limit n — ♦ oo on the 
replicated hamiltonian of RFl.M. The ensuing model is an anisotropic Ising model in 
two dimensions which has short range interactions in one dimension and weak long 
range interactions in the other. The model has been exactly solved. 
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Chapter 1 


Introduction 


The presence of impurities and lattice defects in real materials is an everyday fact of 
life. It is therefore important to assess their relevence and to study the effects of the 
same on the results established in pure cases. The theoretical attempts to study the 
effects of impurities on various phenomena began many years ago but only recently 
has it become possible to develop the necessary theoretical tools to systematically 
analyse the important features. The typical problems include the motion of electron 
in disordered solids, the problem of percolation, spin glasses, etc. [42] 

The presence of impurities have marked effects on the critical phenomena for the 
following reasons. The presence of impurities in materials is equivalent to switching 
on a small perturbation. The response of the system to small perturbations is charac- 
terized by various susceptibilities and correlation functions which become very large 
neaur the critical point. Indeed, they are singular functions of temperatures. Therefore 
even a small amount of impurities can have significant effects. The critical exponents 
which describe the critical bohaviour of various thermodynamic quantities may get 
modified. Even the critical point may disappear, leading to destruction of long range 
order. These effects are nontrivial and are still imperfectly understood. 

We shall mainly be interested in random field phase transitions in which the field 
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conjugate to order parameter is random ( See Ref. ([l]-[3]) for recent review ). The 
random field problem is described by the hamiltonian 

// = - J 53 Si.S, - h,.Si (1.1) 

u * 

where 5, and hi are n-component vectors with 15, | = 1. The random field h, has the 
Gaussian distribution. 

{.srsf) = s‘-%Hi (1-2) 

The case n = 1 corresponds to the random field Ising model (RFIM) while n > 2 
describes the systems with continuous symmetry. 

Ginzburg - Landau version of the model is also very useful in analysis of the critical 
behaviour. In place of discrete spins and discrete lattice sites one has a continuous, 
coarse grained spin variable <^(r) and the hamiltonian is now 

H = jd?? + (Vff - A(r).^(r)| ( 1 . 3 ) 

where 

= Y 

i=l 


vlv$= EDIt)' 


Or=l j=sl 


and the random field distribution is 


(*i(r)A,(r'> = - t‘)S„ . 

» 

The random field model is a good representation of many impure materials. The 
applicability of the model goes well beyond mere ferromagnetic systems. The model 
well describes any solid state system with frozen impurities and point defects which 
locally break the symmetry between ordered states. There are many systems in 
magnetism, ferroclectricity and surface physics where this situation arises [2]. FVom 
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the experimental point of view, Ising antiferromagnet (AF) in an external field He is 
of great significance. Fishman etc. [45] showed that the AF in the presence of external 
field can be mapped on to the RFIM. The AF allows a systematic experimental study. 
The effective random field Hr vanishes with He and can be varied around zero just 
by changing He. This allows one to see how the relevent quantities change with He- 

The field has seen enonnou.s activity in last one and a half decade since the 
pioneering work of Imry and Ma (IM) [4] who pointed out that even small random 
field can have drastic effects on the large scale behavior of the system. Introducing 
basic theoretical concepts, they showed that long range order established in pure cases 
is destroyed by arbitrary small random fields below the lower critical dimension which 
is two for the Ising like systems and four for the systems with continuous symmetry. 

The salient feature of the random field problem is the conflict between the two 
terms in the hamiltonian. The spin coupling J favours allignment of spins while 
random field strength Hr favour.s disoixler. If Hr J, then the system is expected 
to be disordered at low temperatures. However, interesting possibility arises when 
random field is small. 

The usual procedure to determine the lower critical dimension is to consider the 
stability of completely ferromagnetic state to the formation of domains of opposite 
spins. IM proof consisted in showing that the formation of ill-oriented domains is 
always energetically favourable for large enough domain size below the lower critical 
at T = 0. The argument automatically holds for finite T because thermal fluctuations 
favour disorder. 

IM proof is simple and throws light at the heart of the problem, so we shall go 
through it in some detail. Con.sider a bubble (domain) of linear size i in a sea of 
down spins. Let the spins inside the bubble point in the up direction. The domain 
wall energy for Ising like .systems does like as the number of spins that are 
to be flipped is proportional to //"*. For .systems with continuous symmetry one 
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can decrease domain wail energy cost by avoiding a sharp flip at the boundary by 
continuously changing the orientation of the spins along the whole domain. Simple 
arguments show that the energy cost in this case goes like The random field 

energy gain is the sum of L‘^ random variables. The typical veilue of the random field 
energy is therefore ~ with positive or negative sign but it is always possible 

to find a region enclosing an arbitrary point such that Erf > 0. So the total energy 
cost for a bubble of size L is 

E{L) — JL^~^ — for Ising like systems 

= for systems with continuous symmetry 

For large enough L energy cost can always be made negative for d < 2 and d < 4 for 
Ising like systems and systems with continuous symmetry respectively, no matter how 
small the random field is. Tims the ferromagnetic state becomes unstable to domain 
formation for d < 2 for Ising systems and for d < 4for continuous systems. 

Using renormalisation group arguments, IM also verified that the upper critical 
dimension, the dimension above which the exponents are those of mean field theory, is 
six. They suggested e expansion near d = 6 for d = 6—e and carried out the calculation 
to first order in £. Grinstein has [5] considered the 0(£^) terms and found that d is 
replaced by d — 2 in the exponents and hyperscaling relation of the pure system. This 
led to the astonishing idea of dimen.sional reduction whereby critical behaviour of d- 
dimensional random field system is claimed to be same as that of pure system in d— 2 
dimension. .Aharony etc. [6] attempted to prove that dimensional reduction holds to 
all orders in e. The complete calculation has been first done correctly by Young [7]. 
Parisi and Sourl<^[8] avoided the combinatorial difficulty encountered in this approach 
by introducing supersymmetric formulation. This enabled them to present the proof 
in a much more elegant way. 

The dimensionality reduction would imply that lower critical dimension of RFIM 
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is three and that of systems with continuous symmetry is four. For systems with 
continuous symmetry this agrees with IM predictions but is in conflict with their 
results for RFiM. Since the IM predictions are based on handwaving arguments, a 
reconciliation was sought by arguing that IM domain wall arguments give a lower 
bound and there might be a different mechanism which renders ferromagnetic state 
unstable to formation of ill-oriented domains. Attempts were made, [9] and {11], to 
obtain d = 3 as the lower critical dimension and do the e expansion for d = 3+e for R- 
F interface model. It was later pointed out (10) that the qualitative aspect of this work 
has some difficulties. Later work [20]-[21] involving a very plausible R — G procedure 
on R — F interface re-established d/ = 2. The controversy was finally settled when 
rigorous analysis [45] proved that the three dimensional RFIM exhibits long range 
order at zero temperature and small disorder, disproving the dimensionality reduction 
predictions. 

This leads us to question the validity of the analysis which imply dimensionality 
reduction d —* d — 2. The <leri\'ation is based on the assumption that only a selected 
class of diagrams which are most divergent near the assumed transition point are 
relevant in evaluating the critical exponents [6] . These diagrams are called the 
tree diagrams and the approximation amounts to neglecting the thermal fluctuations. 
This is equivalent to solving the local mean-field equations. The difficulty with 
dimensional reduction arguments becomes transparent in Parisi-Sourl^(8] approach. 
They assume that the local mean field equations can be solved and order parameter 
can be expressed as a function of the random field. After this they average over the 
random field configurations. However, they implicitly assume only one solution to the 
local 'mean- field equations. This assumption (as was later pointed by them) is true 
above the transition temperature of pure mean-field problem but breaks down below 
this temperature where multiple solutions are possible. The transition temperature 
in presence of random field lies below the transition temperature of the pure problem. 
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Hence, the implicit assumption which leads to dimensionality reduction breaks down. 
The same will hold in the other approach which in essence is equivalent to Parisi- 
Sourlas approach. 

Schwartz [14] pursued a novel and an interesting idea of constructing an equivalent 
annealed model whose statistical properties are identical to the random field Ising 
model near the critical point. The ensuing model is complicated and has long range 
correlations. But Schwartz claimed that considering the leading behaviour of the 
correlation function, it can be shown that the random field Ising model has the same 
critical behaviour as the pure Ising model in d' = d—2+ri{d). Schwartz indicated that 
this result may be generalised to arbitrary dimensionality of the order parameter. It is \ 
interesting to note that Schwartz reduction gives lower critical dimension in agreement | 
with Imry-Ma result. 

The critical dynamics of the system is aLso drastically modified from dynamics 
of the pure system. E.xpcrirncnts [■2-l]-[28] performed on the diluted Ising antiferro- 
magnets (AF) in uniform field - physical realisation of RFIM - reveal quite unusual 
properties. Since di = 2, AF should exhibit long range order for low temperature 
and small external field. But when cooled in non-zero fields, AF do not exhibit any 
long range order down to very low temperatures even for modest field strength. On 
the other hand long range order established by cooling in zero field persist under 
the application of quite large fields. The history dependent behaviour implies that 
equilibration time is very large. Specifically critical slowing down is so extreme that 
it exceeds the largest observation time. 

The analysis around the upper critical dimension, [30] and [31], indicate that 
the critical dynamics is of conventional type with power law divergence. But this 
is not in accord with experimental findings which show that the releixation times 
are extremely large. To explain the data in terms of power law divergence would 
require unphy.sically large value of dynamical exponent [27], [29]. To account for 
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the extremely large relaxation times Villain [19] and Fisher [18] have suggested the 
activated dynamics where equilibration requires jumps between remote energy wells 
in phase space. The relaxation time, according to Villian and Fisher diverges like 
r ~ exp(^®) where 9 is the exponent that enters the hyperscaling relations. Recent 
experiments [28] have confirmed that the critical dynamics is of activated type but 
the exponent 9 is quite different from the theoretical expectations. 

It is now clear that as far as critical behaviour is concerned, the random field 
fluctuations completely dominate the thermal fluctuations. This has led to the gen- 
eral consensus that critical behaviour is controlled by the zero temperature fixed 
point. The m&in purpose of this thesis is to show the relevence and utility of the 
zero temperature fixed point for both statics and dynamics of the random field sys- 
tems. In the next chapter we introduce the phenomenology of the zero temperature 
fixed point. The fact that there are three independent exponents in contrast to two 
for the conventional systems is stressed. The additional exponent is related to the 
renormalization of the random field. In the third chapter replicated approach is com- 
bined with standard decimation technique to set up renormalisation group flows for 
the random field Ising model (RFIM), which are Migdal-Kadanoff like for dimensions 
d > 2. Considering the RFIM in d = 2 -f c, we find the critical exponents to first 
order in €. 

The systems with continuous symmetry are considered in chapter four. For these 
systems (n > 2), both dimensional reduction and IM domain wall arguments yield 
four as the lower critical dimension. It is thus possible that Parisi-Sourlas reduction 
holds all the way from six down to four dimensions. Indeed, it was claimed in litCTature 
[6] that for n > 2 to leading order in c = d — 4, the critical behaviour of the random 
field model is the same as that of pure model in d = 2-1-e. Fisher [23], however showed 
that these calculations are incorrect. The fixed point found by Young [6] in d = 4 c 
Wcis shown to be unstable to random anisotropies. In chapter four, by combining 



CHAPTER L INTRODUCTION 


5 


decimation and Migdal KadanoiF bond moving approximation we find another fixed 
point in d = 4 + e. Considering the system in d = 4 + c, we find a new set of critical 
exponents to first order in c. 

In chapter five we study the critical dynamics of the RFIM. We show that care- 
ful zmalysis of dynamics at the zero temperature fixed point combines conventional 
dynamics with the activated dyn 2 imics. It is shown that the application of dynamic 
scaling and zero temperature fixed point leads to an everpresent correction to scaling 
which causes the exponent 9 to differ so strongly from the theoretical expectations. 
The last chapter is devoted to one dimensional models in a random magnetic field. 
In section A, we consider the one dimensional Ginzburg-Landau model in a random 
magnetic field. Using the replicated hamiltonian it is possible to reduce the prob- 
lem of calculating the free energy to that of calculating the ground state energy of a 
quantum mechanical problem. In section B, we consider the random field Ising model 
in one dimension. Considering the n-times replicated hamiltonian we concentrate on 
estimating the largest eigenvcilue of the corresponding 2“ x 2" transfer matrix by us- 
ing the variational principle. From the estimated eigenvalue which permits analytical 
continuation, we find the free energy for all temperatures. 



Chapter 2 


Zero Temperature Fixed Point 


2.1 Renormalisation group analysis 

A completely satisfactory renormalisation group treatment of ferromagnets in a ran- 
dom magnetic field is not present although a good understanding of the structure of 

X 

flow equation exist due the pioneering work of Bray and Moore [16] and somewhat 
more general exposition of Berker eind McKay [17]. We shall closely follow Bray and 
Moore and a review article by Nattermann and Rujan [3] in our presentation. 

The model for ferromagnets in a random field is defined by the hamiltonian 

^ = -ill Si ^i-Si ( 2 - 1 ) 

^ (m) ^ *■ 

where the random field has a Gaussian distribution 

= S'‘XH„ ( 2 . 2 ) 

m = HS , (2.3) 

We shall assume an infinitesimal field in direction 1. The indices i, j as usual denote 
the lattice sites and a, whidx refer to the field components run from 1 to n. The 
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cases n = 1, n = 2 and n > 2 correspond to Ising, X-Y and Heisenberg model 
respectively. 

We assume that equilibrium phase transition remains continuous down to T = 0 
for Gaussian random fields. The phase diagram expected for this model is shown in 
figure 1. The existence of three fixed points will be assumed in addition to the trivial, 
high temperature fixed point. 



The three fixed points are 

1. A totally unstable fixed point C aX T = Tg,^ = H = 0, 
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2. A fixed point R&tT = H = 0 and ^ = tn* which is unstable in two but stable 
in in one direction. 

3. A totally stable fixed point O3.tT=^ = H = 0 which corresponds to 
low temperature phase. All flows in the ferromagnetic region terminate at this 
point. 


The critical behaviour of the pure system is controlled by the thermal fixed point C. 
But since random field is a relevent perturbation, flows along the phase boundary 
approach the zero temperature fixed point R. The fixed point R controls the critical 
behaviour in presence of the random field. 

To calculate the critical exponents we linearize the R — G flow close to the fixed 
point R. The eigenvalues and eigenvectors of the linearized R—G flow transformation 
deliver then the critical exponents and the scaling fields. Phenomenological arguments 
concerning the R — G flow suggest t = ^ — w* + j (const) and y as the scaling 
fields. 

On the basis of this picture one can now derive the scaling laws. From dimensional 
considerations the singular part of the fre^as the form 

= (2.4, 

Now imagine carrying out an i? — G coarse graining transformation, with length 
scale factor 6, corresponding to reduction of degrees of freedom by a scale factor 
b^. Introducing exponents yj, and yt which describe the scaling of J, H and t 
respectively, the coarse-grained J', H* and f are given by 

J' = 

t' = tW* (2.5) 





H' = 
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1 

II 

V = Ifpt 

7 = v{2yh -yj-d) 

^ = 2 + d + yj - 

S — xkzM 
d-vh 

y-A-\-d-2yH 


Table 2.1: Scaling Relations 

From equations (2.4) and (2.5) and the fact that partition function is invariant under 
R — G coarse-graining we obtain the following scaling relations for the free energy 
density and the correlation length 

f) = ‘"•'"'/(j*'" f (2.6) 

? ( j. (f (2-7) 

The scaling relations and critical exponents follow in the usual way ( Stanley 1971, [51] 
). For detailed analysis, see Ref.[16]. The relations between yj, pt and ph and critical 
exponents are summarized in the table below 

It is to be noted that there are three independent exponents in contrast to two 
independent exponents in the pure case. The additional exponent appearing in the 
present case is intimately connected to the fact that we are considering the zero 
temperature fixed point ‘R’. In absence of random field the thermal fixed point ‘C’ 
controls the critical behaviour. The analogue of eq.(2.4)is 

F. = Jf{t, j) (2.8) 

where < = j — ( j)* and F, is the singular part of the free energy. At the thermal 
fixed point ‘C’, the exponent yj appearing in eq.(2.6) is zero because J is invariant 
under the R — G coarse graining. But at the fixed ‘R’, it is the ratio ^ which is 
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while the coupling J scales up at the fixed point.From eqs.(2.4) and (2.5) it is clear 
that free energy scales up at the fixed point ‘R’ if yj > 0. This leads to the modified 
hyperscaling law. 

(d — yj)v = 2 — a (2.9) 

The additional exponent rj appearing in Table (2.1) is associated with the scaling 
of the disconnected part of the correlation function 

Giu = = ;3:^Sd,.(r/f) (2.10) 

2.2 Explicit evaluation of critical exponent 

Since the critical behaviour is controlled by the zero temperature fixed point ‘R’ we 
can evaluate exponents simply by working at T = 0. The exponents can be obtained 
from R — G differential equations of the form 


^ = HMw) 

(2.11) 

It 

(2.12) 

f = HJ.M 

(2.13) 


where w = coresponding to the length scale factor b = c*. In general renormaliza- 
tion group transformation generates new terms in the hamiltonian. We shall assume 
that these terms are irrelevent in the R — G sense at least near the lower critical 
dimension. 

The fixed point ‘R’ is obtained from the equation for w. Using eqs. (2.11) and 
(2.12) we get 

^ - fM) {2-14) 

At the fixed point 

dw . _ 

^ 
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giving 


= /2K) 


at the fixed point. Linearizing around the fixed point gives the exponent v, 

yt = ~ = ( 2 - 15 ) 

where the primes denote the derivative of /i(ty) with respect to w. 

The other exponents yj and yh are given by 


Vk = /.(»■) = /2(»') (2.16) 

yn = MW) (2.17) 

Simple domain- wall arguments due to Imry and Ma suggest that the lower critical 
dimension of random field Ising model is two and that of magnets with continuous 
symmetry is four. Near the lower critical dimension the fixed point w* is going to be 
imall. This allows a perturbative treatment of the functions fi{w). 

In what follows we shall use Migdal-Kadanoff bond moving approximation to 
vaJuate critical exponents in d = di + e for the Ising and the Heisenberg magnets in 
random magnetic field. The Migdal-Kadanoff bond moving renormalization group 
“aunsformation was invented by Migdal [39] and reinterpreted using bond moving 
eas by Kadanoff [41]. In Kadanoff interpretation, bonds are moved so that spins 
ty be eliminated by one dimensional decimation. The bond moving approximation 
JK:act for r — » 0 and ^ 0. So we expect that in d = di + e when T = 0 and 

where a > 0; our results are exact. The bond moving predictions have 

'Zy been found to be exact [39] for pure Ising and Heisenberg models near the 
'-f tical dimension. 



Chapter 3 


Random— field Ising Model 

In this chapter we will consider a Random-field Ising Model. The replicated approach 
is combined with the standard spin-decimation technique to set up renormalisation 
group flows which are Migdal-Kadanoff-like for dimensions D greater than or equal 
to 2. 

The partition function of the random field Ising model (RFIM) is 

Z = E “pI? E + ? E Ai'Sil (3.1) 

{S,} ^ (ij) ^ i 

where the random field at the site i has a Gaussian distribution 

{h,hi) = SijHj, (3.2) 

We assume zero uniform external field. We shall use n-times replicated system in one 
dimension which has the partition function (after averaging over the field hi ) ( See 
Appendix 1 ) 

Z" = E«pf?E5r5?;. + #(Esff 

5<^> L *> ^ 

= ( 3 - 3 ) 

c(«) I *,Of a^f3 
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where K =■ ^ and A = ■^. a and are the replica indices and they run from 1 to n. 
The RFIM corresponds to n = 0 in the sense that free energy which is proportional 
to (In Z(k)), can be obtained from /im„_o The idea at this point is to carry 

out the usual renormalization group elimination of degrees of freedom by doing the 
trace over every alternate spin - in this case over the set 3“ - and thus obtziin a flow 
K' = f{K,A,n) and A' = g{K,A,n) and take the limit n — ♦ 0 which gives the flow 
for the RFIM. Indeed, the Taylor expansion about n = 0 is 


IC = f{K, A, n) = K'o + nK[ + o{tP) (3.4) 

A' = Aq +• nAj, + o{rP) (3.5) 

Terms proportional to n and higher order in n on the right hand side of eqs. (3.4) 
and (3.5) contribute o(n^) corrections to the partition function. The free energy is 
given by 


/zfMn— >o( 



and hence is unchanged by these corrections. Therefore the renormalized IC and A' 
are just Kq and Aq respectively. 

We begin by decimating the Ising chain in a random field. The site index i in 
eq.(3.3) is a single integer, and we have spins etc. on various sites. The 

renormalization group procedure would be to do the trace over alternate spins 
S^yS^^Sg^... etc.. If we consider sites 1, 2, and 3, then the spins at site 2 will be 
removed and the partition expressible in terms of 3° and 3^. We thus require the 
relevant bit of Z", which we write as 


^4 


= S exp 
{Sf} 


KY^3t3^ + 3^3^ 

. Of 


exp 


E ^(5fSf + 2S?Sf + s^sl) 


(3.6) 


There are n spins at site 2 which means that we have to sum over 2” configurations. 
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Summing over all these configurations leads to 


ZL 


2 exp 


ji:(s;sf+sssl) 

L 


cosh + X^cosh JflX; - 2<r ' } + .... 

a=l mi a 

_-4(^■-l)(n+l-W)A 


m|m> • • *j| jg* 


,n{n-l)A 


(3.7) 


where 

r=^st+st. 

One can obtain eq.{3.7) by explicitly working out the cases n = 2,3... etc. and then 
generalizing to arbitrary n after observing the pattern that emerges. 

To implement the renormalization- group program, we want to write eq.{3.7) as 


Z; = exp 


Y.K^StSt + E + Sts!) + a 




(3.8) 


To find K% A\ and a as functions of K, A, and n, we evaluate eq.(3.7) and eq.(3.8) 
for three specific situations : 


1. each of {5f} and each of {S^up, 


2. each of {5f } up and each of { 53 } down, 

3. 5P,5f^ up and down, while up and 5?!^ 

The three cases lead to 


gin* +n(«-i)A +a] ^ 2e’*(’*-^>^[cosh(2nife) -I-" cosh{2(« - 2)K}... 

co8h{2{n - 2(iV - l)}i<:] + • • •] 

= 2c’‘t"-')^>l(n) (3.9) 
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2gn(n-I)Ajj . . 

+ • • •] 

2g»(n-l)A 


(3.10) 


g[nA:'+(n-x)(n-4)A'+aj _ 2e('‘-i)(«-4)A [cosh{2(n - 2)K}+ 

{cosh2nJf + (n — l)cosh[2(n — + ... + 

{’‘-'Civ-x cosh[2(n - 2)/'^] Cn .2 cosh[2(n - 2{N - + ... 

= (3.11) 

Considering eqs.(3.9) and (3.10) and dividing one by the other and taking the 
logarithm of both the sides we end up with 

2nIC = In[cosh 2nK +" cosh 2(n -2)K + --- + 

2{n - 2(N - l))K + • • •] 

-In[l + + 

+ • • •] (3.12) 

Eq.(3.12) admits analytical continuation. Expanding R. H. S. of eq. (3.12) in 
power of n we get 

2n/C =s ln[l + n(e^^ cosh 4I( — cosh 8K + 

j. cosh4(n — 1)/^ H ) + o(n^)] 

“ ln[l + n(e^^ - + ••• 

= n{e*^ cosh 4/C — cosh 8/C + • • • + ^ ^ cosh 4(n — 1)K 4 ) 
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- le‘“ + . . . + + .■•) + o(n’) 

Taking the limit n -♦ 0 in eq.(3.13) leads to 

2/r = 2 [cosh(4m/i') - 1] 

msl ^ 


(3.13) 


(3.14) 


It is straightforward to check that for A = 0 eq.{3.14) reduces to the well known 
recursion relation for one dimensional Ising model. 

To find the recursion relation for A, we divide eq.(3) by eq.(3.9) to get 


A{n) 

At this stage we set n = 0 in eq.(3.15) and obtain 

-4 A' _ 4A ^(Q) 

A(0) 

From the expressions for A(n) and C{n) we find 


(3.15) 


(3.16) 


A(0) = 1 


and 


C{0) = ( —1)”*'^^ { cosh 4mK + cosh[4(m — !)/<']} 


in=l 


We may now write eq,(3.16) as 

g4(A'-A) ^ ^ cosh 4m/<' + cosh[4(m - l)iir]} (3.17) 

marl 

Eqs.(3.14) and (3.17) are the renormalization group flows for the one dimensional 
RFIM. 

The series appearing in eqs.(3.14) and (3.17) are divergent and as such are not 
manageable. Using the identity 
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we can recast eq.(3.14) in the form 


K' 


^ r dpe 

J^OO 


j-oo 
1 [°° 
-OO 


-p3 


^ ^^4ntfC’+4mp\/A _j_ g— 4miC+4mp\/A 

n=l ^ 

[ln(l + + ln(l + - 21n(l + 

= /(iir,A) (3.19) 


2 

f dne^^ In 

(1 + ^ g 4 pVS- 4 A'^* 

IIJL 

—OO 

(1 + e.*v^y 


In arriving at the final expression we have used the identity 


n=l ^ 


In an identical fashion, by substituting eq.(3.18) in eq.(3.17) we find 

e4(A'-a) _ ^ 

" 27r«/2 

g4m\/5p-f4(m-l)K' ^ ^4m^/Ap‘~4{m^l)K 


g4(m“l)\/Ap+4mK' ^ ^4(m--l)>/Ap- 

J — OO --.—-I 


AtnK 


+ 

(3.20) 


Now, using the identity 


n=l 1 ® 

we get 

g4{A'-A) ^ r A)e-^dp + r fiiK, A)e-^*dp (3.21) 

J— OO OO 

where 

MK, A) = ^ 

and 

MK.A) = 

We shall work in the regime where T —* 0 and w = ^ is smaU. In terms of K 
and A this implies K -4 0, A -+ 0 such that ^ = tw is small. 
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In the limit of large K and A eq.(3.19) reduces to 

which for small to reduces to 

K* — K{\ — + (exp. small corrections of the type (3.22) 

The random field has the recursion formula given by eq.(3.21). In the limit K —* oo, 
the right hand side of eq.(3.21) is dominated by the contribution from first term of 
fi{K, A) and the second term of / 2 (iir, A). The integration now leads to the recursion 

A' = 2A (3.23) 

which shows the relevance of A and thus at T = 0, the pure fixed point in one 
dimension is unstable. It should be to be noted that, for small u>, the recursion to 
the above recursion involve terms such as which we can safely ignore. 

We now use the Migdal-Kadanoff bond moving approximation to extend the re- 
sults of eq3.(3.22) and (3.23) to higher dimensions. The approximation assumes the 
coupling K appearing in one dimension to be changed to 2^~^K in d dimensions, 
corresponding to accumulated bond strength in a (d — 1) hypersurface - the spins on 
which are removed in a decimation procedure. It is valid at low temperature and small 
random field where an ordered state is assumed so that a simple addition of individual 
bonds can be envisaged. The extent to which the single site interaction A Sa,/? 5f 
should be moved in higher dimensions before decimations is ambiguous. We follow 
Migdal (39] in associating a combined interaction KYlaSiSf + ^f) 

each bond and moving both the spin coupling K and the random field strength A. 
Thus the recursion relations eqs.(3.19) and (3.21) become 


K’ = 


(3.24) 


^ e’ /2(2‘'-'Jr,2‘^-'A)e-*^<fp] 

(3.25) 

where fi and /j are the same functions as before. 

We now demonstrate that, for ty C 1, we get the same structure of flow equations 
as discussed by Bray and Moore [16]. To do so, we consider eq.(3.21) for 2 ui arbitrary 
scale change; i.e.,we consider a factor 6. an arbitrary b would be of order 2" if eq.(3.21) 
is to be generaJized, and then we consider the limit of 5 = 1 + The resulting flow 
has the structure ^ = —AwK where A is an a rbitrary c onst. This holds for d = 1. 
For d > 1 we note that the contribution will come from the factor of atteiched to 
IT, and thus the flow will be of the form 

^^K{d-\-Aw) (3.26) 

As for eq.(3.25), it is straightforward to generalize that, since it leads to A' = 2^A in 
d dimensions and thus 


dA 

dl 


= d A + ( exp. small corrections for tn 1) 


(3.27) 


We have in eqs.(3.26) and (3.27) the form by advocated Bray and Moore [16] with the 
difference that eq.(3.26) is linear in w whereas Bray and Moore argue for a iw* term. 
If we introduce an infinitesimal uniform field H, then it is clear that for T = 0,ty = 0 



(3.28) 


For small tv there will be small exponential corrections. 

In terms of J, Hr and H we can write eqs.(3.26), (3.27) and(3.28) as (omitting 
exponentially small corrections) 


dHn 

dl 


= HrMw) = HR.d/2 


(3.29) 


= Jf 2 {w) = 7(d - 1 - Aw) 


(3.30) 
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dH 

dl 


= Hh{w) = Hd 


(3.31) 


The fixed point is determined from the equation for w = Manipulation of 
eqs.(3.29) and (3.30) leads to 


— = ti;(d — 1 — Aw — d/2) = w{ - ^ ^ — Aw) 


For d = 2 + €, we get 


At the fixed point 


giving w* = 


dl 4 


dw . 

~dl ° 


Linearizing around the fixed point gives the exponent i/ 

y, = l/u = WV{{V>-) - = ^{A) = I 

while the and yj are given by 


(3.32) 


(3.33) 


Vh = = 2 + € 


(3.34) 


yj = /l(u;-) = Mw^) = 1 + €/2 


(3.35) 


to order 0(c). 

Using the results given in the table [2.1], the exponents ti and fj can also be 
computed to 0(c), we find 


T] = 2 + d + yj — 2yh = 1 


^ = 4 + d — 2yj — 2 — c 


CHAFTEJi ^ aAm mfrFlEW miNn MnnRT. 


24 


Schwartz and SofFer [15] have derived the exact inequaility fj t}. To order c this 
is satisfied as equality. Note that ij =. 2r, implies yj = 2 - 7 . This leads us to the 
modified hyperscaling law 

^~oc = d — 2-{-if (3.36) 

to 0 (c). 

Schwartz [14] has argued for exact dimensional reduction for random field magnets 
where the exponents of d-dimensional RFIM are exactly equal to the exponents of 
pure system in dimension <f = d— 2-f ij (d). It is interesting to note that our results are 
in accord with this rule. Indeed ,d d—2+i) implies that d! = l-fc/ 2 . The exponents 
of pure system in dimension 1 -f t/2 are to 0 (c) [52] , 1 /v = e/ 2 , ri = l- c/ 2 , which 
are in agreement with our results. 

Setting d = 4 in eq.(3.32) leads to the flow 

T-'- 

which would give a correlation length ^ ~ e^. Our expression for the correlation 
length differs from the Binder’s [50] estimate of the correlation length ^ ^ 
and that of Grinstein and Ma’s [ 21 ]^ ~ Both these results are based on the 

roughening of the interface due the random field. While Binder’s result is based on the 
discrete version of the Ising model, Grinstein and Ma use the continuum version of the 
Ising model to get their results. Independent calculation of exponents in d = 2 •+• c 
before us has been done by Bray and Moore [16] and T.Nattermann [22]. These 
calculations differ from each other and from us in the renormalistion of J. Bray and 
Moore use first correction to d — 1 in eq.(3-26) as 0(tn*) so as to get in agreement 
with Binder’s correlation length. Essentially same has been done by Nattermann as 
regards Grinstein and Ma’s results. We believe that discrete as well as continuum 
version of the Ising model should belong to the same universality class and hence 
should have the same of critical exponents. Since the calculations based on interface 
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Table 3.1: List of critical exponents in = 2 + e 


o 

II 

- = i 

7^ 

II 

1 

S = oo 

^ = 2 — e 


roughening give different values of critical exponents for the two versions, it seems 
inadequate in explaining the critical behaviour of the random field Ising model. 

In Appendix II, we shall derive the R-G flow using the conventional technique 
where one renormalises the unaveraged hamiltonian and keeps track of the probability 
distribution of the random field and the spin coupling. 




Chapter 4 


Random-field Heisenberg Model 


For magnets with continuous symmetry (n > 2) both dimensional reduction and 

Imry - Ma domain wall arguments yield four as the lower critical dimension. It is 

thus tempting to speculate that Parisi - Sourla’s dimensional reduction holds all the 

way from six dimension down to four dimension. It was claimed in literature [6] that 

for n > 2 to leading order in « = d — 4 the critical behaviour of the random field 

model is the same as that of pure model in d = 2 + c. Fisher [23], however showed 

that these calculations are incorrect. The fixed point found by Young [6] in d = 4 + e 

was shown to be unstable to perturbation by random anisotropies. These authors 

start with n - times replicated hamiltonian of the non-linear a model with random 

•«* — ♦ 

m^netic field. The replicated hamiltonian is then expanded about x = = 1 

when ail the replicas are aligned along the same direction where a and ^ are replica 
indices. We believe that it is not consistent to expand about x = = 1, when 

one is interested in the perturbative fixed point in Indeed, for ^ = 0 the replicas 
are independent and {Sa-S^)oi ^0 = 0. For ^ small, (5o,.5/j)a//3 is uot near one but is 
near zero. It is perhaps this fact which disables one to find the fixed point in d = 4-l-e. 

The X — I seems to be the ground state configuration but it is to be noted that 
the random field related contribution to the partition function from this configuration 
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is of O(n^) which has no physical significance. Consequently, our aim will be to avoid 
the use of replicas in our derivation of the flow equations. 

In our analysis below, we avoid replicas and use simple decimation and Migdal - 
Kadanoff bond moving approximation to study the Heisenberg model near d = 4 + e 
The partition function of the Heisenberg model is 



where the random magnetic field has a Gaussian distribution. 


(if 

(if) = HSu. 


(4.2) 

(4.3) 


We shall assume an infinitesimal field in direction 1. 

We begin with af = 1, the Heisenberg chain in a random magnetic field. The 
renormalisation procedure would be to trace over alternate spins i.e.Sf , 5'!', ... etc. at 
various sites. If we consider the site 2* — l,2t, and 2* 4- 1; the spins at site 2i will be 
removed and partition function is expressed in terms of and After this 

coarse>gratning we replace by Sf. We thus require the relevent bit of Z which 
we write as 


Z' = f d52,exp(~(^&*_i.52i + 52i.52f+i) 

JSj, I 

+ 2h2i.S2i + A2t+l-^2t+l)] 


2^n/3 


4»/a-i (f + ■52f+i -f ^l) 
(f |^2»+l + + ^l) 


n/2-1 


I 


exp{^(ft3i-1..92,_l + A2i+l-'S2t+l)] 


(4.4) 


Obviously tho hauiiltoniau does not preserve its forni after coarse-graining. However, 
we shall work in the regiirw where T 0 and ti> = ^ is very small. In this limit the 
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average angle between two neighbouring spins Si and Si+i is very small and we madce 
the usuaJ spin wave expansion, 

Si.St+i = cos0,,,+i = 1 — (higher order terms in . 

In the limit of T — » 0, we can approximate eq.(4.4) by 


Z' w exp 


S^i^i + Sj.+i + 


*2i 


+ const i exp ^{hi -l.S 2 t-l + ^2»+1->5^2«+1 ) (4.5) 


We are working in the regime where where we may expand 4 821-1 + S 2 i+i + ^ as 


follows (treating b — the small expansion parameter ) 


J 

T 


§ 21-1 + 5j,+i + 


/J2, 


J 


; I 2 + 252»-l.‘52t+l + 2-^.(5'2 ,-x + ‘52.+i) + 


^2».^2» 
J2 


1/2 


.h 2 i 


h 2 i.h 2 i 


1/2 


y 4 - ^2t_1.2i+l + 2-^.02i-l + <?2.+l) + 


7* r \ 


2J ( . ^2»-l,2«+l , ^ /Pt , rt \ , ^2i-b’2i 


r 


4 +27-(^-> + ^“«>+ 4J2 


2/ /. ,, h2i.h2i. 1 T* 


) 77^2«.(^2i-l + ‘52 «+i)X 


8P ' 4J 


^ ” g^2i-l,2i+l{^ 

•«* 

h • h ’ 

+ higher orders in b + terms of the sort 6 ^ x 

8 




(4.6) 


where x, y, and z are integers. New terms in the hamiltonian are generated. We 
assume that they are irrelevent in the R—G sense at least near the lower the critical 
dimension. Therefore we may write eq.(4.6) as 
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£ 

r 


5*21-1 + ^a.+j + 


>2i 


7 <'Lu.«{i-y^} + 


42 


~.(5ai-i + S 2 ri+i){l - + irrelcvent terms 

2 T 


h 2 i.h 2 i ■ 




8J2 

i* 

^•{^ai-i + S 2 ,+i){l - - ^— ■ } + irr. terms (4.7) 
By substituting eq.(4.7) in eq.(4.5), we end up with 

h2i-«2ii^ « h2i-h2f^h2i 


2; 


exp 


J 

2 T 


{I - + {1 - ^}^.(5a.-a + 52m) 


8J2 


8J2 ^ 2 T 


1 - - ■ 
+ ^(^21-1*52,-1 + hai+j.^ai+i) 


(4.8) 


Similaxiy we may decimate the spin at site 2i - I to arrive at 


Z' = exp 


J ,, h2,_2*h2,_2 , ^ ^ , f, h2<-2*^2»-2i ^2 j- 2 /cf • C \ 

/b2i-3*'i>2i-J + (1 gj2 ;-2^-(‘^2i-3 + *2i-lj 


1 


+ ^(^21-3*52,-3 + ^2,-1 *52,-1) 


(4.9) 


After replacing ^ 2 i~i,^ 2 i -3 and by §i,§i-i and §i+i respectively, we get the 
reRormali2!ed J' and as follows 

= (4.10) 

Af = At-. + Y (‘ - + % (' - + 0(6^) (4.11) 

After taking the average over Af’s in eqs. (4.10 ) and (4.11) by using eqs.(4.2) and 
(4.3), we end up with 


J 

H'^ 2 H 


+ 0 (w^) 


~ 2 V 8J2 y 

(,_M)+0(»^) 


(4.12) 

(4.13) 
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£ 

T 


Sii-i + + 




fi ^2t-h2is 

- ^jTf + 

1* ■*+ -4 

+ 52i+i){l - + irrelevent terms 

|j2^^ 

^•{52i-i + 52 ,+i){1 - + *»'»*• terms (4.7) 


By substituting eq.(4.7) in eq.(4.5), we end up with 


ZL 


exp 


1 o ? , <1 ^2i-^2ii^2» 

2r' 8J2 


■}52i-i-52i4.i + {1 


gjj }2^.(52,-i + Sa+,) 


1 


+ ^(^2«-I'S'2i-i + ^ 2 <+ 1 - 52 i +1 


(4.8) 


Similarly we may decimate the spin at site 2i -- I to arrive at 


ZL 


exp 


J f. ^3M-3-ti3i-3 ^ S . fi ^2i-2-^2i-2i^2«-2 ^ \ 

8«/^ /^2i-3‘^2»-i "i" SJ^ ^ 2T* ' 


i 


+ .jj,(^2»-3-^'2i~3 + Aai-i-Sai-i) 


(4.9) 


After replacing 1 ,^ 2, -3 and .^ai+i by and <^+1 respectively, we get the 

renormalized J* and ft! as follows 


fi - ?km) 
V 82^ ; 


(4.10) 


+ 0(6’) 

/ 

‘f = *?.-■ + Y (‘ - ) + % (1 - + 0(6’) (4.11) 

After taking the average over hf*s in eqs. (4.10 ) and (4.11) by using eqs.(4.2) and 


(4.3), we end up with 


‘'"21 8J2 


^ + 0(w* 


) 


(4.12) 




(n + 2}ffl 

16J2 


j + 0(w*] 


(4.13) 
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in 4 - 2)H% 
12J» 



I 'nfortunatfly, r<*n«rmiUisation Ifads to a new term 


(4.14) 


V - ikfhrl,) = iM 




in -f 2)ffi 


4^2 


+ 0(u?^} 


which was zero to start with. The ratio 


m = 


{hfhf) 


increa.'«*s after each R-G iteration. A proper choice of this ratio would make it 
i nvariant under the renormalisation group transformation at the assumed fixed point 
tv* s= w = u'*. This is done by taking 

(l - + 0(t»*) (4.15) 

i n the original hainiltonian. After the renormaiisation, using e<|.(4.10) and (4.11) we 

*' = f 1 - + 0(w*) (4.16) 

Hi = W'kf) = fl - + 0(w‘) (4.17) 

"The ratio 

<1 8 j ' 

remains the same after R-^G iteration, (ofcourse, w also changes after each R-G 
deration, but if we are able to locate a fixed point in tn, then v/ = w = w*). The 
?aq.{4.14) is however modified to eq.(4.l7). 

These relations hold in one dimension for T -» 0 and small w. Using Migdal - 
fleetdanoif bond moving approximation in which we replace J 
%dl H -* there is an ambiguity in moving single spin interactions, however, 
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we can use the same procetlure as described in the last chapter) the recursion relations 
are mwlified to 


1 

ft 


(4.18) 

«' = 2'‘H f 1 


(4.19) 

«i = 2!w„( 


(4.20) 


If i« is small the bond moving should be a very good approximation because for most 
of the ras<»9 neighbouring spins will be aligned in the same direction and therefore 
bond moving costs small energy. 

An arbitrary scale change b is given by li = 2^ after N R - G iterations.In the 
limit & = l+Sl wc get the structure of flow equations as discussed by Bray and Moore 
(16] up to 0(«;*) 

^ (5 - + 0 („<) = HkMw) (4.21) 

~ ss J ^cf ~ 2 - nAu/^ + O(te^) = J /^(w) (4.22) 

^ 4- 0(w*) ^ H fsiw) (4.23) 

where A is an arbitrary constant. Henceforth we shall drop terms involving higher 
order corrections in w. 

The fixed point R is determined from the equation for tn = ^. Manipulation of 
eq8.(4,21) wrd (4.22) yield the recursion 

~ =r m(/,(ttf) ~ Mw)) = w 4- 1^— ^ . (4.24) 

For d 5s 4, eq.(4.24) leads to the correlation len^h ^ ~ exp Putting d = 4 4- e, 
we get 


dw / e , (n-2)Aaj’\ 
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/\t thr |K»int 

du' . 

l»M**f* ~ 0 

giving 

ty* ss 7"' *ii j ^ ■ 

Linearizing arnunrl the fixf*d point yields the exponent y, = 

V. = i = 

= »!••{- (n + 2) Aw* + 2nAw*) = {n - 2)Aw*^ = e (4.25) 

while exponents yjr and y* are given by 

,. = /,K)=4 + <-^< = 4 + «[l-^^l (4.26) 

yj = MW) = MW) = 2 + «/2 - = 2 + |(1 - ^1 ( 4 . 27 ) 

From cq».(4.25)to (4.27), we can easily find ail other critical exponents. The values 
of critical exponents are listed in the table (4.1). 

Schwartz and Soffer (15) have derived the exact inequality q < 2iy . To order c this 
is satisfied as equality.This would also imply that d-dependent hyperscaling law is 
modified to *2 ~ or = i^{d ~ 2 + iji)’ in presence of the random magnetic field. However 
there is no exact dimensional reduction, rf -♦ d~2-f-9(= 1 + for tbe exponents, 

since for d = 1 4* the exponents of the pure problem to 0(e) are, 1/u — ^e 
and ri ss both being in disagreement with the result given in table (4.1). For 

i» ~ 2 the fixed point remains undetermined. This means that special treatment is 
required for this case as Is the case for the pure model without the random field. 







Chapter 5 


Dynamics 


In this chapter we turn our attention to dynamics of the random field Ising model. 
The critical dynamics of the RFIM is drastically modified from the dynamics of the 
pure system. Experiments performed on the diluted Ising antiferromagnets (AF) in 
uniform field - the physical realisation of the RFIM - reveal quite unusual properties. 
Since the lower critical dimension of the random field Ising model is two , AF should 
exhibit long range order for low temperatures and small external fields in three di- 
mensions. But when cooled in non-zero fields, AF do not exhibit any long range order 
down to very low temperatures even for small field strengths. On the other hand, long 
ramge order established by cooling in zero fields persists under the application of quite 
large externsd fields. The history dependent behaviour implies that it takes very long 
time for the system to come to equilibrium. Particularly, the critical slowing down is 
so extreme that it exceeds the largest observation time. 

The analysis around the upper critical dimension (da = 6) indicate that the criti- 
cal dynamics is of the conventional type with power law divergence [30]. To explain 
the data in terms of the power law divergence would require very large value of the 
dynamical exponent [27]-[29]. This value of the dynamical exponent seems unphysi- 
cal. To account for extremely large times, Villain [19] and Fisher [18] suggested the 
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activated dynamics where equilibration requires jumps between remote energy wells 
in phase space. 

The dynamics of the model has been the subject of extensive experimental in- 
vestigation [24] - [27] in an effort to understand whether the critical dynamics is 
conventional [30]-[31] or of the activated variety [18]-[19]. The initial experiments 
were not able to resolve the difference between the two forncis of dynamics ( the span 
of frequency and equivalently correlation length was not large enough) - activated dy- 
namics could not be differentiated from the conventional dynamics with a very large 
dynamical scaling exponent(of the order of 10). In the more recent determination of 
a.c. susceptibility Nash, King and Jaccarino [28] covered a wide range of frequency 
and could definitely assert that the dynamics was of the activated variety, although 
the relevent exponent found by them did not agree with the theoretical predictbns.In 
this chapter we apply dynamic scaling and zero temperature fixed point idea to the 
RFIM dynamics somewhat more carefully than has been done before and find a cor- 
rection to scaling which serves to reconcile various apparently contradictory features. 
We show that the exponent 0 can be directly determined from the imaginary part of 
the susceptibility. Rirther it is demonstrated that this 0, which is larger than that 
quoted by Nash et al from an analysis of the real part of the susceptibility, is in 
agreement with there result once the “correction to scaling ” has been taken in to 
account. 

The two approaches to critical slowing down have been 

1. conventional dynamic scaling [11] where the relaxation time r ~ ^ being 

the correlation length. The exponent z is found by constructing a diagram- 
matic perturbation theory in which the temperature induced fluctuations are 
completely ignored. 

2. dynamic scaling with activated dynamics [19]-[18] being the basis for critical 
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slowing down, giving r ~ e 2 rp(^*), where 6 is the scaling exponent for the 
lemf^r&iure field. In this approach the random field is important for producing 
the melastable configurations , while temperature is what causes the hopping 
between metastable states, giving rise to the activated dynamics. 

We show that a very straightforward application of the renormalisation group 
transformation fat the relaxation time with the dynamics governed by the zero tem- 
perature fixed combines the two pictures given above. We recall briefliy the salient 
features of the zero temperature fixed point. In this case the hamiltonian is no longer 
scale invariant and under the scale transformation by a factor b , is transformed from 
H to where 9 is an exponent associated with dimensional reduction. Alter- 
natively the temperature can be considered to have transformed and we have the 
transformation law 

r = b'^T (5.1) 

For the conventional fixed point of the critical dynamics without the random field, 
sa 0 and the temperature is unchanged. 

Now turning to the dynamics we note that under a scale transformation b, the 

correlation leni^h goes from ( to ^/b and the frequency scale goes from a> to b~*u), 

where z is a dynamic scaling exponent and hence for the characteristic time t 

T(?,r) = 6'T(|,n = >‘’-(f.3’r*) (5.2) 

Setting 6 = 

Hi, T) = VHl, TC^) = r /(^) (5-3) 

Jf ^ 5 = 0 , it is immediately clear that conventional dynamic scaling is obtained. 
For « 0 , the functional form of /(T/C®) can be crucial. The two factors in eq.(5.3) 

i* and f{TJiH cmrespond to the ojnventional dynamics and activated dynamics 

respectively. 
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We need to determine the function /(T/^®) . To do this we note that this part of 
T is clearly produced by the fluctuations induced by the temperature. The random 
field induced fluctuations in r are contained in the part In calculating f{TlC^) 
, the role of random magnetic field is to produce several metastable states and the 
temperature provides the noise that leads to hopping between the minima. The 
effective temperature is and is very small for large This dynamics can be 
taken to be governed at temperature T by a Fokker Planck equation of the form 

^ = V.{PVF) + kTV^P (5.4) 

at 

where P is the probability of a particular field configuration at a time t and the scalar 
field ^(z) is the coarse-grained order parzuneter for the RFIM. The potential 
is the one determining the deterministic drive in the Langevin Equation for ^(z). 
The potential V{4>) is characterised by the existence of several metastable minima. 
In zero dimensions ^z) is an algebraic variable, the Fokker Planck equation is one 
dimensional and it is a well known result that the largest relaxation time in the system 
has the structure ezp(AV/r) where A V is the typical height of the potential barrier. 
In non-zero spatial dimension, the Fokker Planck equation is infinite dimensional 
corresponding to infinite number of Fourier components of ^(z). The relaxation rate 
is still of the form exp(AF/r), provided tunnelling path exists [32]-[33] connecting 
the different roetMtable minima. The typical potentied barrier is of order imity at 
these low temperatures and hence at temperatures of Tf(^ 

leading to 

T = exp(Ti(^)') (5-6) 

where To and Ti are system dependent constants. If the metastable minima are ig- 
noredy which is equivalent to taking d = 0 ( temperature does not scale ), the relax- 
ation rate would be t = with z obtained from a perturbative expansion around 




D = 6, leading to 

^ = 2 + 2j7 (5.7) 

Our principle result is that given in eq.(5.6). We first point out that if r is taken to 
be dominated purely by activated dynamics, then eq.(5.6) implies that 

t 

T = TQexp{T^{^Y^^f) (5.8) 

where 



For reasonable correlation lengths it is clear from eq.(5.9) that < 9. It should be 
noted that the coefficient z/r^ of the correction term in eq.(5.9) is not necessarily a 
small term and hence the correction to the exponent 9 can continue to be effective 
well into the critical regime. In the experiments of King et al [26], the largeness of 
z/ti leads to the conclusion that the critical behaviour may be governed by activated 
dynamics or by conventional dynamics with a faurly large value of z. The experiments 
of Nash et al [28] are carried out in a more asymptotic range and hence the critical 
behaviour of activated dynamics is seen. However, the exponent 9 is lowered from its 
zero temperature fixed point value. 

If the data were analysed in terms of conventional dynamic scaling on the other 
hand, eq.(5.6)wo«id imply 

T = T„(|)'-" (5.10) 

where 

j.,, = 2r + 9T.(i)' (5.11) 

If Ti is a number considerably smaller than unity, this would be reasonable for small 
but would yield large values of z*// for large correlation lengths. The typical z 
obtained by King et al [26] was around 14. We believe that the same situation holds 
for the experiments of Geshwind et al [29]. 
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Wc now a^ldress the cjuestion of applying dynamical scaling when the relaxation 
time is given by eq.(5.6). To set up dynamic scaling we need to express in terms 
of the relaxation time t. Noting that exp^® will dominate numerically over the 
critical region, we immediately find 


For the temperature and frequency dependent susceptibility x ( or any response 
function that is measured ), the critical behaviour will be given by 




In 

/ 


(5.13) 


The exponent i determines the zero frequency, that is the thermodynamic behaviour. 
If we now pass to the extreme non-thermodynamic Umit i. e. ^ oo and the frequency 
alone determines the response, we have for a finite response 


xKT = r.) ~ (/ii|w(|lnw|)S|)' 


(5.14) 


For X 0 (which is what happens for this response) 


x(w,!r = Tc) In (in jmfnu?{)»j) 

A frequency dependent response aiT — Tc will be a complex quantity with the real 
and imaginary part following the Kramers-Kronig relations. This is due to causality, 
which can be ensured by replacing by -iu in eq.(5.14). Hence at T = Tc, 

Xi + *X2 = Xo (fn(-ta;[/n(-ta;)]t)) ^ (5.15) 

If we approximate z/0 ~ 1 and consider a frequency range where \Itm\ is considerably 
larger than x/2, then the ima^nary part xa of x has the form ( for x 0) 

y. 1 — - (5.16) 

In |a;/nu;] 
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Clearly m order of magnitude smaller than the real part xi, which scales as 
In 1 . The prediction of eq.(5.16) should be easily verifiable experimen- 

tally. 

We will now demonstrate that in the hydrodynamic regime { small ^ ), the rise of 
Xs with increasing ^ can be quantitatively obtained from d 3 mamic seeding considera- 
tions. The hydrodynamic re^me in the present case is obtained for (f/ln ^lnuj\ C 1. 
Writing the scaling form of x as in eq.(5.13) with w replaced by — ia;, we find that 
the leading terra in x In the hydrodynamic regime is 

ci t» i6+x tx+B 

y ^ fX > _ tX > _ > I *” ^ (K 17') 

ln{—iu>lnu>) ^ fn(wlnu?) — i| ln{u;lnu)) 2 (/n|a;/nu?|)* 


leading to 

Y2 -j - 

{In |u;lru<;|) 

For ar — ► 0, we have the prediction 


(5.18) 


X3 ~ 


M§ 1 


(5.19) 


At a fixed frequency, the imaginary part of the susceptibility gives a direct determi- 
nation of 6 and using the data of Nash et al., we find 9 = 1.25. 

We now show that this 0 and that obtained from eq.(5.9) are consistent. To do 
80 we note that the range of wo/w for the experiment of Nash et al is 10* to 10*. The 
middle of the range is at about 10® * which leads to of ~ 10 from eq.(5.12). 

The earlier analysis of King et al, which was in terms of a power law yielded an 
exponent of about 14. FVom eq(5.11), this corresponds to ti w 0.8 . If we now use 
tbi * value of T\ in eq.(5.9) , we com® up with 9 — 9cjf ^ 0.20. With the 9eff as found 
by Nash et al as 1.05 ± 0.16, we find 9 « 1.25 in agreement with that found from the 
imaginary part the susceptibility. 

We end by providing the results for the case where the wavenumber is finite. The 
finite frequenty, finite wavenumber response function x(6‘*'>^) should be amenable 
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to neutron scattering experiments. The dynamic scaling requirement would give the 
form 

= (5.20) 

Clearly for 0, /(fcO must tend to unity. The more interesting limit is clearly as 
-+ oo i. e. approaching the critical point = oo) at finite k and w, In this limit 
clearly 

= (5.21) 

which implies that the k dependent relaxation rate at ^ = oo goes as 

Q{k) ~ fc*exp(-A:"®) 


This relation can be probed by neutron scattering at finite k and w as the critical 
temperature is approached. 



Chapter 6 


One Dimensional Models 


6.1 Ginzbtirg-Landau Model: Equivalent Quan- 
turn Mechanics 

While the random field Ising model has seen a lot of activity, the low dimensional 
Ginzburg-Landau model has not attracted much attention. In low dimensions (i.e. 
!> = I or 0 ), there is no phase transition (not even at T = 0 as there is in the 
absence of the random field ) and any calculation of a thermodynamic quantity, e.g. 
the free energy, must produce a result which is valid over a large temperature range. 
For a one dimensional Ginzburg-Landau model without the random field, this has 
been done by Ferrell and Scalpino {34] and by Balian and Toulouse [35] who showed 
that the problem of calculating the free energy of the model could be reduced to the 
problem of calculating the ground state energy of a quantum mechanical problem. 
For the zero dimensional Ginzburg-Landau model in a random external field , Bray 
and McKaae [36] used the replica trick to evaluate the partition function. In this 
section, we treat the the one dimensional Ginzburg-Landau model in a random ex- 
ternal field by a combination of the two techniques. We use the replica trick and 
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then Jin&lytif ally continue to imaginao-y values of the spatial coordinate to set up the 
equivalent quantum problem. The ground state energy now leads to the free energy. 
The free energy functional for the scalar Ginzburg-Landau in an external field is 


ft 

I dx 
H 


+ i 


dM{x) 




dx 


+ ^M*{x) — M{x)h{x) 
4 


( 6 . 1 ) 


The random fieldl h{x) has the probability distribution 

h^{x)dx 


( 6 . 2 ) 


(6.3) 


P(h) = exp|-/ 
and we are interested in the partition function 

Z{h{z)l = j£)[M]exp-|;|-jf F{x)d3 
A change of scale 4 = Mffi on the field variable yields 

Z[&(*)] jDl^jexp ^^j^L[<l>]dx\ 

Witrll 

L14] = |ro^*{®) + ^ ~ 

The differential element i?[^] is defined as usual as a limit of a product over discrete 
values of x and o is some lattice spacing. Rom the partition function the thermody- 
namic free energy can be obtained by avera^ng over h and is gven as 


(6.4) 


(6.5) 


G = {lnZ)h 


( 6 . 6 ) 


Differentiation now gives the specific heat. The difficulty of calculating Z{h{x)\ and 
then performing the h average is circumvented by the replica trick, whereby we study 
the n-times repBcated system 


Z' = j n D{<l>„\ex^-L] 


(6.7) 
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where 


i=k5:«>(x)+i5: 

jp|t jQp 


d^(x)Y , 1 „ _ 

-_J 4 . _ ^ ^ (s) _ _ ^(x)A(x) ( 6 . 8 ) 

The replica index a runs from 1 to n. Now the h average can be performed to yield 

in = j U 1'' dx\lroj:^Ux)+ 


Ot I 


a^a(x) 


dx 


+ T53^iW-EW^)*W 

^ iTK 


exp 


/ A®(i)dx 
J " 4A 


f w f 

/lle^P -I 

OrSSI i 


where 


4# «% 


1 

2 


d<f>a(x) 


dx 


e// 


+fE«W-| 




12 


(6.9) 


( 6 . 10 ) 


If we analytically continue the coordinate to imaginary valuey. e. consider x -+ it, 


then Ltff becomes the lagrangian for a quantum mechanical problem with 


i = i 

T 

or 

# 

\d4c,] 

. J 

1 

_ _ 

E 

or 

’d<l>a 
dt . 


I'-o x; + J E w - 

4 ^ ^ Of 


V(^a) 


( 6 . 11 ) 


We can consider the n variables as the coordinates of a particle moving in a 
n-dimensional potential well The quantum hamiltonian associated with the 

lagrangian is given by 


H 


1 V- ^ ^ 


iroE^J(x) + fE^lW-f 


12‘ 




( 6 . 12 ) 


Denoting the eigenvalues of H by £«, the partition function {Z ) at temperature T 
in the limit L-* 'oois given by 


(Z") = 


(6.13) 
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where €o(n) is the ground state energy of the hamiltonian. 

One can easily calculate the thermodynamic free energy density as 

-y(/n2) =riim^^-|^lnr (6.14) 

Ju n-*0 dfi 2(2 


We now proceed to find the ground state energy of the hamiltonian of eq.(6.12) 
by writing it as 


If - .15]^ + 
2^d4>l^ 




A 

T 




L or 


( 6 . 15 ) 


Our strategy is to split H as 

H=-'£H^ + H, ( 6 . 16 ) 

Or 

where 

rr 1 . 1 ^2/ ^ /ft t7\ 

~ + 2^0<f>ai^) + 4 ( • ) 

and 

= (6-18) 

C(0 

The hamiltonian /.' j* will now be repl 2 w:ed by equivalent harmonic oscillator hamil- 
toniaa //o«, for the ground state, such that 


= + ( 6 - 19 ) 

where ^f>o is present only if r© < 0 

One can now diagonalize Ea ^oa+Hi to find the approximate ground state energy 
for H. Clearly, the task would be to find as accurate an answer for Wo, as possible. 
Accordingly we have to worry about two distinct regimes (i)ro > 0 and (ti)ro < 0. In 
the case when Tq > 0 we write 


Ifoa 




( 6 . 20 ) 
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where 

= '•0 + ( 6 . 21 ) 

an(i{<>* ) is the expectation value of in the ground state energy of the hamiltonian 
//oa- The self consistency is assured by treating the part ^ ^ 

perturbation and requiring that the first contribution to the ground state energy 
vanish. This leads to the condition 




( 6 . 22 ) 


which is automatically satisfied by our chosen hamiltonian Hoa, and hence from 
eq.(6.2I), u;| is given by 


w ‘ = To 4- 


ZuT 

2t4?a 


(6.23) 


Clearly, to* is independent of the index a. For r > 0, there is only one positive root 
of eq.(6.23) in W(,. We immediately see that eq.(6.15)can now be diagonalized as 


2 " 2 V ' 


(6.24) 


where w* * w* - 2n| for a = I (say) and equals w* for the remaining (n - 1) values 
of a. We can now immediately write down the ground state energy as 


co(n) = - 2n— ] + -(n - l)u; (6.25) 

where w is the solution of eq.(6.23). The free energy follows from eqs.(6.14)and (6.25) 


BS 


p ^ It[^ - - ^InT 

2 ^ rJ 2a 


(6.26) 


Tq <.Q i In this case there are several extrema and — 0 corresponds now to a 
maximum. To find the minima we differentiate V{<j>) with respect to (j>a. 


vm 






L « 


(6.27) 



nfAPTK H 6. O SEMMF.SSIDSAL hfODELS 


47 


liT A ^ 

V'li) = roO, + (6-28) 

■* ^ 0=0 

Setting V"(0) = 0 and considering the replica symmetric solution we get 

- ro<i> * - 2n^^* 


which implies 


/— ro + 2nAfT 

^==^V — — • 


Using the shifted variable 0 = - 4>*, we may write the hamiltonian H as 

= 


1 1 , 6?iA 

”2 ? ^ "F 

’a ^ »/J 


'£,^l + \luT{-ro + 2n^) + 




(6.29) 


where 


and 


Wo, = - + \/uT(-r-M + 


^ Ot0 


(6.31) 


We have dropped those terms in Hqc which give 0(n*) contribution to the ground 
state energy. By re-adjusting the terms in Hoa » we may write it as 


E 


Oa 


2 dtpl 2 2 




— k (6-32) 

4ur ^ 

We now replace Hoa by the equivalent harmonic oscillator 5o« for the ground state. 
We identify 


^ + i(2|r„| + 'fcMd’H + ^ 


(6.33) 
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and calculate ail the average values with respect to the ground state wave function 
of the hainiltoriian. Accordingly 


where 


and 


^ulrolT , 

(6.34) 

UJa 

(6.35) 

(0“) = -3C,^^!^(1+«/3) 

(6.36) 

= ^(l+25-f5V3) 

(6.37) 

S = 

'*'Of 

(6.38) 

nT 

= 21rol + 

(6.39) 


The correction factor S is small in all regimes of uT and hence the consistency con- 
dition of 

(^a) = 

and 


( 6 . 40 ) 


(’/’D = 

sets C, = Cl a 3 . We determine (d’») = ('/'*) W) « independent ot a) from 

eqs.(6.35), (6.38) and (6.39) as 

a 1 , 9ur[rol(0^))^ 

^ " 2llro| + (3ur/2)(W L ^ I2hl + (3uT/2)((i/>W^\ 

and as a byproduct we write down the approximate ground state energy of ifoa as 


( 6 . 42 ) 


If, , 3uT,, 3}-'- y i^oityy u /q_ 

€b * ^ [2|ro! + J - 22ro + Zl2uT{tk^) iuT' 


9 uTlroKC^^))^ rl 


( 6 , 43 ) 
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by using e<is.(6.33) and (6.34). This will correspond to the case when there is no 
random magntftic field. 

An immetliate question is, how good is this approximation ? For uT = 1, we 
can compare with the numerical results of Ferrell and Scalpino [34]. The agreement 
is excellent, in particular, the position of zero of c(ro) (the most sensitive test of 
accuracy of the approximation ) is reproduced for within 5%. 

Having determined the equivalent harmonic oscillator, we now turn to tackle Hj. 
We are faced with the problem of finding the gfbund state energy of 

H = (6-M) 

a ^ ap 

This can be achieved by the orthogonal transformation 

ijf = 0 t 


where O is the orthogonal matrix 



With this transformation 



The ground state energy is now easily seen to be 



uT\ro\m^ nrl 
- 2nA 4ur 


(6.45) 


(6.46) 


(6.47) 
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from which (tee energy foiiows aa 



(6.48) 

_ TV A ^ ri T 

2 2uf 4ti 2a 

(6.49) 

where 



(6.50) 

and is giwen by eq.(6.42). 



Discussion: We note that as compared to the no extemed held case of Ferrell and 
Scalpino, free energy is decreased for all values of tq. For ro > 0 the decrement 
decreases as ro is increased. If we define 


AF(ro) = F(ro)-F(ro)U.>« 

then AF(ro) is negative definite, goes to zero as jroj -* oo and has a peaJc at negative 
value of tq. It is instructive to define the entropy function S’(ro) = j from 

which we see that because of the effect of random field , the entropy is increased 
for al! values of ro. We conclude with a discussion of specific heat whidi is given by 
C ss — In the absence of random field, there is a peak in the specific heat at a 
negative value of ro- In the presence of random field the peak still exists. In fact, 
the peak is more pronounred as A-terra makes a positive contribution. For positive 
values of r© and r© > 1, the specific heat goes to zero as it must reproduce the A = 0 
limit. The behaviour of the specific heat in the one dimensional random field should 
be accessible experimentally and therein lies the utility of the above results. 



6.2 Random Field Ising Model 


The random field Ising modtd in one dimension has attracted considerable attention 
in recent years [46] - [49]. Though there is no phase transition in one dimension, 
the presence of random field introduces new features which are absent in the pure 
problem. In this section we will estimate the free energy of RFIM in a nonperterbative 
manner by using the replica trick. 

The hamiltonian of the random field Ising model in one dimension is 

H = -J'£S,Sn.,-'£hiS, (6.51) 

i i 

where h, is random with Gaussian probability distribution 

{kiki) = SiiH% (6.52) 


. We axe interested in obtaining the thermodynamic free enery 

(FH = {lnZ), (6.53) 

As mentioned in the previous section, this may be obtained by using the replica trick, 
whereby we study the n-times replicated system. The average partition function of 
the n-tim^ replicated system is given by 

(Z"> = E“pf?E5fSf« + #E5“5?] 

{Sf} V"* / 

ss y] exp f/(:V5f5?„ + A5:SfSf| (6.54) 

(S,.) \ ie »/! / 

I JNf^ 

where K =: f and A = - 3 # . 

Using the familiar transfer matrix method, the evaluation of the partition function 
( 2 ®) of one dimensional chain can be reduced to the computation of the largest 
eigenvalue of the 2® x 2® matrix T(s,sO. The transfer matrix T(s,/) has the form 

T(s, s') = exp{~ Z S^S^] explK £ exp[| Z 

A I I I 
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where S = ,5") and S' = (S'^S'*, ,5'") 

In the therniodynajnic limit, the free energy is given by 

F = r^l„„ (6.56) 

The procedure relies on the fact that the form of A(n) would permit analytical contin- 
uation in n. Since we are unable to diagonalise the matrix T{s, s') for any arbitrary n, 
we calculate the approximate largest eigenvalue by using a method which is equivalent 
to using trial wave function in quantum mechanics. 

We consider the trial column matrix 


gn»a/2 \ 

g{n-2)*a/2 


{"Cielements) 

g(n-4)*a/2 

(’‘C'aelements) 

g(n-6)^A/2 

• / 


(6.57) 


This choice of elements of the column matrix is motivated by the fact that it is 
the eigenvector of r{a,s') which corresponds to the largest eigenvalue A(n) in two 
extreme limits, namely when 


1. K ss any arbitrary value and A = 0 


2. K = 0 and A « any arbitrary value. 

The largest eigenvalue may now be computed as 

X(n) sas (in ihe Dirads bra - ket notation) 

(flfa) 


(6.58) 
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We have iiunuTScally checked this for n = 2 and n = 3 when it is possible to diagonalise 
the transfer nratrix T . For all T and A that were numerically checked, we found 

A(2) - M 

— < 0.18. 

By taking the expectation value of T with respect to ja), we end up with 
{«!^la) = E £ ’^Cmr”'CrCo Ci_x'"(7ic"^ + . . . 

I lms:0 

^n-mj^^mQg4ia-)g{n-W-2m)K-g(n-2m)*A] 

^{n-2l)\/iShi+{n-^)K j riQ^ fn-mQmQ^^ 

/ Yl 1^0 

Ct-rC\e*^ + . . .""”‘Co"*Cje^^^') 


n-m 


/ Y' 

J sfisft y 


^{n-20N/Sfc|+{n-2l)ic|^n^^^^^ 


(6.59) 


where 


/tmn “ ("“"CrCo + "‘*^Ci_rC,e"^' + 

To evaluate the series eq.(6.59) we consider the function 

/««(x) « (1 + (6.60) 


It is easily verified that coefficient of in the Taylor expansion of /mn is Ann- 
Therefore the the coefficient in the sum ^Cmfmn is equal to '£,^ "Cm An».But 

III3S0 mss0 

ifisO ^ 

^ / -JK-2VS». (6.61) 
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Therefore the roelficieiil of x* is 


™ ^ g2ir-2\/AAiy 

Substituting this in Ktj.(b.59), we get 

(airja) = / ^ g2v^M')’‘-' 

/ V jsiQ ^ 


^^tsTKhi ^ g2/v*-2V^fciyj 
- ^n^^^g-/f-2>/Sfc,+2VAAj ^g/f+2v^Ai+2>/AA2^’*-' 

-2%/2rha gK-2\/Afcj-2N/Sjia^^ 

fdhxdh^ 


f * ** 

J \/ 5 r\/i 




glir+'2\/S(fcj+fcj) g-*/f+2V*A(fci — fcj) 

■w '■ V ■■■ 

^g-K--2i/S(A,-/ii) ^ gK--3^{h,+fc,)l“ 


(6.62) 


In a similar way we obtain 


(o|a) = e"*^+"CiC<"-’J*'^ + + + , 

J y/w 


(6.63) 


We thus find 

J ^1 4- g-y’f2VS(A3-Ai) ^ g-/f-2>/A(A3-/>i) ^ gy-2v^(Aa-l-fti))» 


J ^f,-h^(^g2y/Eh ^ e-2'/aA)n 
By using the eq(6.56) and eq(6.64), we obtain the approximate free energy as 

P <j> f dfii dA2 ^_A;-ftg |jj/g/f+2s/a(Aj+a2) g-/if+2\/a(A2-Ai) ^ 

“ J JxJv 


(6.64) 


_|,g-lf-3\/3i(&a"“Ai) ^ gK—iy/S(h7+hi)'^ 
J y 5r 


(6.65) 
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It b instruftivr to runsid.^r the free energy and entropy in the limit of T = 0. The 
frn* energy is fount! out to bi* 



while entropy is zero at T = 0. This is in disagreement with other published 
results which give a finite value of entropy at T = 0. But it is to be noted that we are 
taking the distribution of the random field to be continuous. The non-zero value of 
entropy in the discrete model conies from the fact that the random field exsists over 
a set of measure zero, which is not true for the continuous distribution that we are 
working with. 

One can improve the approximation by introducing a variational parameter a and 
then maximizing A{n,a) with respect to a. This can be done by modifing the trial 
column matrix a. We take the trial column matrix as 

^ gn*ora/2 ^ 
g(n-2)*aA/2 


("Cielements) 

g(n-4)®Qa/2 


("Cjelements) 

g{n-6)®aA/2 


(6.67) 


The largest eigenvalue can now be estimated in a similar fashion as done before 


in this section. We find 

A(n,a) = - 


( 6 . 68 ) 





M 


where 

= e^‘+2v^I^(fc,+A,)^g-A'+2V2|iA(fc,-fci) 

^ g-K-2y'^^p^(fca-fci) g/f-2y''S^(Aj+ii) (6.69) 


and 


B( Ji , a) = (6.70) 

The value a = a that maximises A(n,a) in the limit of n 0 is given by 

/ = / ~e-'‘'^lnB(A,a) (6.71) 

J s/t ^/n da ' J y/n da 

The free energy is now easily seen to be 

F = T I In A{K, A, a) 

J ylT 

^fi-'‘’lnS(A,a) (6.72) 

At T =s 0, the expression for free energy is simplified to 


I e-‘V/i) -2^ (6.73) 

where 6 = ® determined from 


(6.74) 



Appendix I 

Jk. 


Thie Replica Trick 


In addi tion to the therinal averaging one must confront the task of averaging over the 
random impurities in the disordered systems {For details see [53]). Basically there 
axe two types of random impurities which are classified as annealed impurities or 
quenched impurities according to the way they are distributed in the host system. 

Let 5, ( I =s 1 to N) denote the statical variable (spins in magnetism ) and {x} 
represent the impurity configtiration. In annealed system the impurities are essentially 
in thermal equilibrium with the host system and so are averaged in the same manner 
as the statistical averages; for instance, the free energy of the system is 

F « ■-kTHZ{x]U (AI-1) 

where 

Zix} ^ E expl-Ii{x,Si}fkT] (AI.2) 

is.) 

In case of quendied impurities , the impurities in the host system are fixed while 
the statistical variable fluctuate. Based on the argument due to Brout [54], one can 
deflne a meaningful average whereby one finds the free energy or some other extensive 
variable for a given impurity configuration and then average over it. For instance the 
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average free energy is given by 

F=-ibr{InZ{x,5}]<,„ 

The fact that In Z must be averaged rather than Z makes the problem of quendxed 
impurities much more intractable than its annealed counterpart. To average over 
In Z the unaveraged partition function has to be evaluated. Technically it is nearly 
impossible as one has to specify the value of typically 10^4 random variables for a 
macroscopic system. One way to circumvent this difficulty is to use the replica trick 
which is based on the identity 

lim ^ = Inx (AI.3) 

For positive integer n, one writes the partition function of n identical replicas of the 
system. 

Z’(x) = n ZM = n E «pi-iH(j;,sf)i = £ exp[-uTi;ff(^.sf)i 

oral a»I 

where Za(x} is the partition function of the replica. 

At this stage, one averages over Z’^(x} and expresses partition function in terms of 
effective hamiltonian ffe// which contdns no disorder and is translationdly invariant. 
The disorder in the system is mMiifested as some sort of an interaction between the 
replicas. The partition function may thus be written as 

(z-> = y; «p[-«^,(n)] (Au) 

W) 

To be specific we consider the RFIM, where the hamiltonian is 

(«i) 

Tbe partition functicna of n-repHcas is 

z'{fc.)=exp(iE(iEsr5i'+E‘.5r) (a“) 

^ (Oral 
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Averaging over the random variables h[s, which we assume to have a Gaussian dis- 
tribution leads to 


{ 


fZSfSf 
r w) 


S 

T2 






(ALT) 


where //#* is the random field strength. 

Assuming that the analytical continuation to non-integral values of n is possible, 
one can evaluate the free energy by using the identity given in eq(AI.3). Therefore 


the free energy is 




from which all other thermodynamic quantities may be evaluated. 


(AI.8) 



Appendix II 


Renormalisation Group For RFIM 


In chapter 3, we derived renormalization group flow for the RFIM in one dimension 
by using the replica trick. In this appendix, we show that one gets the same structure 
of equations if one avoids the replicas and uses the conventional techniques where 
one renormaliscs the unaveraged Hamiltonian and keeps tracks of the probability 
distribution of the random held and spin coupling J. 

The partition function of the random field Ising model (RFIM) is 

Z = E «p f ? i: SiSi + i E ( AII.1) 

V • / 

where the random ma^aietic field hw a Gaussian distribution. 

(AII.2) 

{hi) = H (AII.3) 

We shall assume to be infinitesimally small. As has been done in Chapters 3 and 4, 
we begin with ~ 1, the Ising chain in a ramdom magnetic field. The renormalisation 
procedure would be to trace over alternate spins, . . . etc., at various sites. 

If we consider the site 2i — l,2i, and 2i + 1, the spins at site 2i will be removed and 
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partition function will he expressed in terms of S^i-i and Sj.+i. After this coarse- 
graining we replace by 5,. We thus require the relevent bit of Z which we write 

m 


^ ~~ 5-f + — (+h2,-52,)] (AII.4) 

“ ^xp{~{52,_i + Saf+i) + -^1 + exp[-^{52i_i + 52, + 1 ) - —] (AII.5) 

~ exp{— {52 ,_i52, + 52ii52i+i) + ^ {52i_i + 52,-4.i + a)] (AII.6) 

The three utiknowns in eq.( AILS) can be evaluated by equating eq.(AII.5) and eq.(AIL6) 
for three specific cases. 

1- = ^ 3—1 ~ 1 


2. 52 , 4 .} “ ‘^3,-1 = ““I 

3. 52,4 j = — 52,-1 = I 




After a simple algebraic mimipuIatioQ we find 

We are parliculariy interested in the re^pme where T 0 and w 
this limit we jfiad ( apart from exponential small corrections in u? ) 


A 


(AII.7) 

(AII.8) 
^ < L In 


Jl = J-jIM 

*' = 5*- 



Siiiice we axe eliminating spins at even sites, the field at odd sites is not affected. It is 
clear that the field at the odd sites will also contribute to the renormalisation of the 
random field. After some thought one gets the following expression for J‘ and W. 


Jl = J- (AII.9) 

K = *K-. + jAs + (An.io) 

Eqs.(AIL9) and (AII.IO) yield ( after taking average over fii’s) 

= J - -^Hr + 0(e-'/'^) (AILll) 

v2t 

«K = |«K + 0(e-'''^) (An.l2) 

+ (AIL13) 

Unfortunately, renormalisation leads to correlations between neighbouring spins. In- 
deed 

4 

which was zero to start with. 

The ratio 

{hihi) 

increaaa alter each R-G iteratioil. A propar choice of this ratio would make it m- 
variaat under the reoormalisatiou group transformation at the assumed fixed pomt 
» = «.•. This i. done by taking (hih,+.) = iBl (we shaU drop 0(e-*/'^) corrections 



M 


heriraft«*r ) in the original Hamiltonian. After the renormalisation, using eqs.(AIL2) 
and {Ail. 10) wc get 

X' = = -Hfi (AIL14) 

//« = {hrM) = 2Hl (AIL15) 

The ratio 

(W "4 

remains the same after R-G iteration (of course, w also changes after eadi R-G 
iteration, however, if we are able to locate a fixed point in w, then «;' = «? = w*). 
The eq,(AII.12) is however modified to eq.(An.l5). 

These relations are true in one dimension for T -+ 0 and small w. Using Migdal - 
Kadanoif bond moving approximation in which we replace J J, Hji -* 
and H ( there is an ambiguity in moving single spin interactions, as such 

we use the same procedure as described in the chapter 3 ) the recursion relations are 
modified to 

J' = 2‘*-V(i--^) (AIL16) 

H' r: 2'^H (AIL17) 

ff'fi = 2^^j, (AIL18) 

If to is smdll the bond moving should be an excellent approximation because most of 
spins will be aligned in the same direction and bond moving costs small energy. 

An arbitrary scale dhange b is ipven by 6 =; 2^ after N R — G iterations.In the 
limit h =s I + we get the same structure of flow equations as discussed in chapter 3. 

(AIL19) 

dl ^2 

^*J(d~l-A«o) (AII.20) 

m 
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(AII.21) 


where A is im ajrhitr&ry constant. 



Appendix III 


Anisotropic Ising spin system in 
two dimensions 


In this appendix we study a two dimensional anisotropic spin system on a rectangular 
lattice of N rows and M columns with spins 5f , * = 1, and a = 1, 2, ....N located 
at the intersection of i*^ row and column. The hamiltonian of the system is 

= f;sf5f-AEEsf (AIII.1) 

IWlCRXl |‘=sl a=5l 


where J is the nearest neighbour coupling for along the rows while ^ is the weak 
long range coupling for spins along the columns. Though this model has nothing to 
do with the random fields we include it here because it has the same form as N times 
replicated hamiltonian of the random field and is amenable to exact treatment. 

The partition function of the system is 




isn 


M ff 


M N 


M N 


+ 5^ sfSf(9fEE5? 

* tael flrsel tasl orrsl 


The model can be solved by using the important identity 


(AIII.2) 
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exp 




2NT 






T ‘ 


The partition function may now be written as 

W) 


/ V Ir \ ^ r I M N M T \2 M N 

^ = (^) /,n".E«p ?EEs?si.-ivE#-EE( 

\4XI / /-so , i-ic-i ,_j , =10,-1 


JiAi + h 


S? 


T 

(AIIL3) 

It is clear thafijA,’ is the effective magnetic field at the column. For a uni- 
form system X is expected to be a constant. Though this assumption is suggestive 
and intuitively obvious, we are unable to make it more rigorous. This assumption 
considerably simplifies our task. First of all, we evaluate the quantity 


f / ^ ^ ^ 7,A + A 

2'= E«p ^EEs?sr„+EE^sr 

{5J*} 1*1 «»k| 1*1 Or*l 


(AIIL4) 


This can be done by using the transfer matrix method. It is to be noted that Z* is 
nothing but the partition function of the N times replicated Ising model in a uniform 
field. Therefore 

= (Ani-5) 

where pi and pa tl** eigenvalues of the Ising model in a uniform magnetic field. 
W Pi > Pa, we can neglect the contribution of pa as compared to pi in the Umit of 
M oo . So, finally we end up with 


where 


Pt 


\2HT)Lo. 




(AIII.6) 

(Ani.7) 





m. 


md 

4W = ^-Tln(h 

la the limit of M — » oo , iV -♦ oo; the integral can easily be evaluated by the method 
of st«*p<^t deacent. One easily finds that the free energy density is 

/ = Ail) (MILS) 


vffhere A is the value that minimizes A(X). The minimum condition 


dA 

dX 


lx=X = 0 


implies 


sinh(^* 

8inh^( ^ * |^ " " ) + e”"^] 


(Ain.9) 


The value of A for which A( A) is minimized has the physical meaning of magnetization 
density and hence we write eq. (AIIW) as 


- - (AIII.10) 

"" (8inh*(=^) + c-^]V2 

The solution of eq.( AIH.IO) is most easily found by the standard graphical technique. 
There axe three solutions for T<T,md only one for T > Z. The critical temperature 


is given by the equation 

ll 

Cb 

1 

(AIII.il) 

Clearly m must be small for h - 0 and T near Tc for T < Tc- 

In this repme 

eq.(AIII.iO) simplifies to 



m 


(AIII.12) 



(AIII.13) 

Or 


(AIII.14) 
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which giv«>s the value of the exponent 0 = 1/2. One can also easily find all other 
exponents. Thetw are found to have the mean-field values. The effect of the short 
range interaction i« completely suppressed in the critical regime though the critical 
temperature is found be very sensitive to the changes in J as can be easily checked 
from eq.(AIIMl). 
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